home *** CD-ROM | disk | FTP | other *** search
/* clvrstat.c v.1.10 This program parses a CLOVER channel status file. It prompts for the name of the HAL channel status file to be parsed. Normally this has the format nnnnnnnn.tst, where nnnnnnnn is the input file's time tag. Output to the console and the file "nnnnnnnn.sum" is various summary statistics calculated from the status file. Output to files called nnnnnnnn.MY and nnnnnnnn.HIS separates the input data into space-delimited MY and HIS form suitable for plotting or further analysis. The program also makes a total of seven MY and HIS gnuplot setup files. You can ignore these if you don't use gnuplot. See the November 97 QEX paper for the names of these setup files and a description of their content. Preface all comment lines after the first in the HAL input file (as many as you want) with a "#" (pound) symbol. Note that many comments in the source code use the C++ comment-line marker ("//"), which is accepted by most modern C compilers. If your compiler doesn't, enclose comments in the usual C comment symbols. Note also that the program doesn't have all the crash tests of a commercial application. Drive at your own risk! Plot-formatting file names are now DOS-compatible: "mysnrtpt.plt" "hssnrtpt.plt" "myphsecc.plt" "hsphsecc.plt" "mydfreq.plt" "hsdfreq.plt" "myhspwr.plt" Written by K. Wickwire. 5 May 1996. Modifications. v1.01 09.05.96 Output to input.MY and input.HIS. Make plotset files. v1.02 17.05.96 Read TO and FROM calls from input file. v1.03 01.06.96 Launch gnuplot on exit with (non-ANSI) exec() function. (Not in this version.) v1.04 08.06.96 Re-organize plots: my/his snr&tput, phs&ecc, dfreq; make corresponding plot setup files. v1.05 03.01.97 Log time in input data as HH/MM/SS for gnuplot 3.6beta. v1.06 13.04.97 Re-label "TO" and "FR" as "MY" and "HS". v1.07 19.04.97 Extend treatment to FEC status data (only one callsign in file). v1.08 14.08.97 Clean up for public distribution. Use 8.3 plot setup names. v1.09 16.08.97 Set rmargin 5 characters wide ("set rmargin 5"). v1.10 09.09.97 Modify to handle DSP4100 channel stats (auto power, etc.). */ #include <stdio.h> #include <string.h> //For string-processing. #include <stdlib.h> //For atoi(), etc. #include <math.h> //For sqrt(). #include <time.h> //For date-time stamp. #define numchars 70 //Maximum input-file line length. int main(void) { short len,i,MY_call_found,HS_call_found; short max_MY_tput,max_HS_tput,max_MY_SNR,max_HS_SNR; short no_MY_lines,no_HS_lines,no_MY_XX,no_HS_XX; short no_MY_num_ecc,no_HS_num_ecc; //Nos. of numerical (non-XX) ecc percentages. short no_MY_BPSM,no_MY_QPSM,no_MY_8PSM,no_MY_8P2A,no_MY_16P4A,no_MY_16PSM; short no_HS_BPSM,no_HS_QPSM,no_HS_8PSM,no_HS_8P2A,no_HS_16P4A,no_HS_16PSM; long sum_MY_bias,sum_HS_bias,sum_MY_tput,sum_HS_tput,sum_MY_SNR,sum_HS_SNR; long sum_MY_dfrq,sum_HS_dfrq,sum_MY_phs,sum_HS_phs,sum_MY_ecc,sum_HS_ecc; long sum_MY_pwr,sum_HS_pwr; long sum_sq_MY_bias,sum_sq_HS_bias,sum_sq_MY_tput,sum_sq_HS_tput,sum_sq_MY_SNR; long sum_sq_HS_SNR,sum_sq_MY_dfrq,sum_sq_HS_dfrq,sum_sq_MY_phs,sum_sq_HS_phs; long sum_sq_MY_ecc,sum_sq_HS_ecc; long sum_sq_MY_pwr,sum_sq_HS_pwr; double avg_MY_bias,avg_HS_bias,avg_MY_tput,avg_HS_tput,avg_MY_SNR,avg_HS_SNR; double avg_MY_dfrq,avg_HS_dfrq,avg_MY_phs,avg_HS_phs,avg_MY_ecc,avg_HS_ecc; double avg_MY_pwr,avg_HS_pwr; double sd_MY_bias,sd_HS_bias,sd_MY_tput,sd_HS_tput,sd_MY_SNR,sd_HS_SNR; double sd_MY_dfrq,sd_HS_dfrq,sd_MY_phs,sd_HS_phs,sd_MY_ecc,sd_HS_ecc; double sd_MY_pwr,sd_HS_pwr; time_t now; //For date-time stamping of files. struct tm *date; char s[80]; int getline(char *line, int max, FILE *iop); //Prototype for line-getter. char fname0[20],fname1[30],fname2[30],fname3[30]; char *ext1 = ".sum",*ext2 = ".MY",*ext3 = ".HIS"; char MY[10]; //MYcall ("MY"). Handles AB1CDE-12. char HS[10]; //HIScall ("HS"). char temp_hr[3],temp_min[3],temp_sec[3]; //For time conversion. char line[numchars]; char firststr[30],secondstr[10],thirdstr[10],fourthstr[10],fifthstr[10]; //Data-line strings. char sixthstr[10],seventhstr[10],eighthstr[10],ninthstr[10],tenthstr[10]; FILE *input,*output1,*output2,*output3,*output4,*output5,*output6,*output7,*output8,*output9,*output10; printf("Enter name of CLOVER input file: "); //Get name of input file. scanf("%s", fname0); if((input=fopen(fname0,"r")) == NULL) //Open input file for read only. {printf("Cannot open %s\n",fname0); return 0;} for(i=0;i<30;i++) //Initialize file names. {fname1[i] = 0; fname2[i] = 0; fname3[i] = 0;} strncpy(fname1,fname0,8); //Copy 1st 8 chars. of input file name into fname1. strcat(fname1,ext1); //Add ".sum" to fname1 to make output file name. strncpy(fname2,fname0,8); strcat(fname2,ext2); //Make ("nnnnnnnn.MY") output file name. strncpy(fname3,fname0,8); strcat(fname3,ext3); //Make ("nnnnnnnn.HIS") output file name. if((output1=fopen(fname1,"w"))==NULL) /*Open .sum file for statistics.*/ {printf("Cannot open main output file."); return 0;} if((output2=fopen(fname2,"w"))==NULL) /*Open .MY file for graph data.*/ {printf("Cannot open .MY file."); return 0;} if((output3=fopen(fname3,"w"))==NULL) /*Open .HIS file for graph data.*/ {printf("Cannot open .HIS file."); return 0;} if((output4=fopen("mysnrtpt.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open mysnrtpt.plt file."); return 0;} if((output5=fopen("hssnrtpt.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open hssnrtpt.plt file."); return 0;} if((output6=fopen("myphsecc.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open myphsecc.plt file."); return 0;} if((output7=fopen("hsphsecc.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open hsphsecc.plt file."); return 0;} if((output8=fopen("mydfreq.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open mydfreq.plt file."); return 0;} if((output9=fopen("hsdfreq.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open hsdfreq.plt file."); return 0;} if((output10=fopen("myhspwr.plt","w"))==NULL) /*Open plot file for parsed data.*/ {printf("Cannot open myhspwr.plt file."); return 0;} len = 0; //Initializations. MY_call_found = HS_call_found = 0; no_MY_lines = no_HS_lines = 0; no_MY_XX = no_HS_XX = 0; no_MY_BPSM = no_MY_QPSM = no_MY_8PSM = no_MY_8P2A = no_MY_16P4A = no_MY_16PSM = 0; no_HS_BPSM = no_HS_QPSM = no_HS_8PSM = no_HS_8P2A = no_HS_16P4A = no_HS_16PSM = 0; sum_MY_bias = sum_HS_bias = 0; sum_MY_tput = sum_HS_tput = 0; sum_MY_SNR = sum_HS_SNR = 0; sum_MY_dfrq = sum_HS_dfrq = 0; sum_MY_phs = sum_HS_phs = 0; sum_MY_ecc = sum_HS_ecc = 0; sum_MY_pwr = sum_HS_pwr = 0; sum_sq_MY_bias = sum_sq_HS_bias = 0; sum_sq_MY_tput = sum_sq_HS_tput = 0; sum_sq_MY_SNR = sum_sq_HS_SNR = 0; sum_sq_MY_dfrq = sum_sq_HS_dfrq = 0; sum_sq_MY_phs = sum_sq_HS_phs = 0; sum_sq_MY_ecc = sum_sq_HS_ecc = 0; sum_sq_MY_pwr = sum_sq_HS_pwr = 0; max_MY_tput = max_HS_tput = 0; max_MY_SNR = max_HS_SNR = 0; time(&now); //Get & write clock time. date = localtime(&now); strftime(s,80,"%d.%m.%y %H:%M:%S",date); //Get MY and HS calls. while( (len != -1) && ((MY_call_found == 0) || (HS_call_found == 0)) ) { for(i=0;i<numchars;i++) //Initialize line[]. {line[i] = 0;} for(i=0;i<10;i++) //Initialize strings to be parsed. {firststr[i] = 0;secondstr[i] = 0;} len = getline(line, numchars, input); //Get numchar characters from line. if(len < 39) continue; //Skip first label line. MAKE SURE IT HAS < 39 CHARS! if(line[0] == '#') continue; //Skip line if it starts with "#" (comment line). for(i=0;i<numchars;i++) //Replace commas with spaces in data. {if(line[i] == ',') line[i] = ' ';} if(MY_call_found == 0) { sscanf(&line[0],"%s %s",&firststr,&MY); //Get MY call. MY_call_found = 1; continue; } if(HS_call_found == 0) { sscanf(&line[0],"%s %s",&firststr,&HS); //Get HS call. HS_call_found = 1; } } //End (first) while. rewind(input); //Go back to beginning of input file. printf("Statistical Summary of CLOVER channel status data:\n\n"); fprintf(output1,"Statistical Summary of CLOVER status data.\n"); fprintf(output1,"Input file is [%s].\n",&fname0); fprintf(output1,"MYcall (\"MY\") = %s; HIScall (\"HS\") = %s.\n",MY,HS); fprintf(output1,"Summary date: %s\n\n",s); //Time stamp results. fprintf(output2,"#Status data for %s.\n",MY); fprintf(output2,"#Input file is [%s].\n",&fname0); fprintf(output2,"#Summary date: %s\n",s); fprintf(output2,"#(ecc = \"XX\" written as \"110\")\n\n",s); fprintf(output2,"#H/M/S mod b tp SNR df phs ecc pwr\n"); //Column labels. fprintf(output3,"#Status data for %s.\n",HS); fprintf(output3,"#Input file is [%s].\n",&fname0); fprintf(output3,"#Summary date: %s\n",s); fprintf(output3,"#(ecc = \"XX\" written as \"110\")\n\n",s); fprintf(output3,"#H/M/S mod b tp SNR df phs ecc pwr\n"); //Column labels. while(len != -1) //While loop through data lines of input file. { for(i=0;i<numchars;i++) //Initialize line[]. {line[i] = 0;} for(i=0;i<10;i++) //Initialize strings to be parsed. {firststr[i] = 0;secondstr[i] = 0;thirdstr[i] = 0;fourthstr[i] = 0; fifthstr[i] = 0;sixthstr[i] = 0;seventhstr[i] = 0;eighthstr[i] = 0; ninthstr[i] = 0;tenthstr[i] = 0;} len = getline(line, numchars, input); //Get numchar characters from line. if(len < 39) continue; //Skip first label line (assumed < 39 chars). if(line[0] == '#') continue; //Skip line if it starts with "#". for(i=0;i<numchars;i++) //Replace commas with spaces. {if(line[i] == ',') line[i] = ' ';} sscanf(&line[0],"%s %s %s %s %s %s %s %s %s %s", &firststr,&secondstr,&thirdstr,&fourthstr,&fifthstr, &sixthstr,&seventhstr,&eighthstr,&ninthstr,&tenthstr); //Line format: //time ,call ,mod ,b,tp,SNR, frq,phs,ecc,pwr (b=bias; 1=robust,2=normal,3=fast) //203808,KB1JY ,BPSM ,2,02,031,-3 ,024,000,000 (tp=thruput; 0-99 bytes/s) //203808,W1IMM ,BPSM ,2,02,028, 003,031,050,000 (mod=BPSM,QPSM,8PSM) if(strcmp(secondstr,MY) == 0) //Parse MY-lines. { no_MY_lines++; if(strcmp(thirdstr,"BPSM") == 0) no_MY_BPSM++; //Count modulation choices. if(strcmp(thirdstr,"QPSM") == 0) no_MY_QPSM++; if(strcmp(thirdstr,"8PSM") == 0) no_MY_8PSM++; if(strcmp(thirdstr,"8P2A") == 0) no_MY_8P2A++; if(strcmp(thirdstr,"16P4A") == 0) no_MY_16P4A++; if(strcmp(thirdstr,"16PSM") == 0) no_MY_16PSM++; sum_MY_bias += atoi(fourthstr); sum_MY_tput += atoi(fifthstr); sum_MY_SNR += atoi(sixthstr); sum_MY_dfrq += atoi(seventhstr); sum_MY_phs += atoi(eighthstr); sum_MY_pwr += atoi(tenthstr); if(strcmp(ninthstr,"XX") == 0) { no_MY_XX++; //Count # times ECC capacity is exceeded. if(strcmp(MY,HS) == 0) no_HS_XX = no_MY_XX; //To make FEC counts correct. //Replace "XX" with "110" (arbitrary replacement). ninthstr[0] = '1'; ninthstr[1] = '1'; ninthstr[2] = '0'; } else //Update ecc percentages. { sum_MY_ecc += atoi(ninthstr); sum_sq_MY_ecc += atoi(ninthstr)*atoi(ninthstr); if(strcmp(MY,HS) == 0) //File has FEC stats. { sum_HS_ecc += atoi(ninthstr); sum_sq_HS_ecc += atoi(ninthstr)*atoi(ninthstr); } } sum_sq_MY_bias += atoi(fourthstr)*atoi(fourthstr); sum_sq_MY_tput += atoi(fifthstr)*atoi(fifthstr); sum_sq_MY_SNR += atoi(sixthstr)*atoi(sixthstr); sum_sq_MY_dfrq += atoi(seventhstr)*atoi(seventhstr); sum_sq_MY_phs += atoi(eighthstr)*atoi(eighthstr); sum_sq_MY_pwr += atoi(tenthstr)*atoi(tenthstr); if(atoi(fifthstr) > max_MY_tput) max_MY_tput = atoi(fifthstr); //Update max tput. if(atoi(sixthstr) > max_MY_SNR) max_MY_SNR = atoi(sixthstr); //Update max SNR. sprintf(temp_hr,"%c%c",firststr[0],firststr[1]); //Get time tag. sprintf(temp_min,"%c%c",firststr[2],firststr[3]); sprintf(temp_sec,"%c%c",firststr[4],firststr[5]); //Print time-tagged data to "MYcall" file. fprintf(output2,"%s/%s/%s %s %s %s %s %s %s %s %s\n", temp_hr,temp_min,temp_sec,thirdstr,fourthstr,fifthstr, sixthstr,seventhstr,eighthstr,ninthstr,tenthstr); } //End MY-line processing. if(strcmp(secondstr,HS) == 0) //Parse HS-lines. { no_HS_lines++; //Count modulation choices. if(strcmp(thirdstr,"BPSM") == 0) no_HS_BPSM++; if(strcmp(thirdstr,"QPSM") == 0) no_HS_QPSM++; if(strcmp(thirdstr,"8PSM") == 0) no_HS_8PSM++; if(strcmp(thirdstr,"8P2A") == 0) no_HS_8P2A++; if(strcmp(thirdstr,"16P4A") == 0) no_HS_16P4A++; if(strcmp(thirdstr,"16PSM") == 0) no_HS_16PSM++; sum_HS_bias += atoi(fourthstr); sum_HS_tput += atoi(fifthstr); sum_HS_SNR += atoi(sixthstr); sum_HS_dfrq += atoi(seventhstr); sum_HS_phs += atoi(eighthstr); sum_HS_pwr += atoi(tenthstr); if(strcmp(MY,HS) != 0) //Different MY and HS calls. { if(strcmp(ninthstr,"XX") == 0) { no_HS_XX++; //Count # times ECC capacity is exceeded. //NOTE: ALREADY SET ABOVE TO 110 IN FEC DATA! ninthstr[0] = '1'; ninthstr[1] = '1'; ninthstr[2] = '0'; } else //Update ecc percentages. { sum_HS_ecc += atoi(ninthstr); sum_sq_HS_ecc += atoi(ninthstr)*atoi(ninthstr); } } //End if(strcmp(MY,HS) != 0). sum_sq_HS_bias += atoi(fourthstr)*atoi(fourthstr); sum_sq_HS_tput += atoi(fifthstr)*atoi(fifthstr); sum_sq_HS_SNR += atoi(sixthstr)*atoi(sixthstr); sum_sq_HS_dfrq += atoi(seventhstr)*atoi(seventhstr); sum_sq_HS_phs += atoi(eighthstr)*atoi(eighthstr); sum_sq_HS_pwr += atoi(tenthstr)*atoi(tenthstr); if(atoi(fifthstr) > max_HS_tput) max_HS_tput = atoi(fifthstr); if(atoi(sixthstr) > max_HS_SNR) max_HS_SNR = atoi(sixthstr); sprintf(temp_hr,"%c%c",firststr[0],firststr[1]); sprintf(temp_min,"%c%c",firststr[2],firststr[3]); sprintf(temp_sec,"%c%c",firststr[4],firststr[5]); //Print data to "HIScall" file. fprintf(output3,"%s/%s/%s %s %s %s %s %s %s %s %s\n", temp_hr,temp_min,temp_sec,thirdstr,fourthstr,fifthstr,\ sixthstr,seventhstr,eighthstr,ninthstr,tenthstr); } //End HS-line processing. } //End while-loop through lines. if((no_MY_lines < 2) || (no_HS_lines < 2)) {printf("Not enough data for statistics.\n"); return 0;} avg_MY_bias = (double)sum_MY_bias/no_MY_lines; //Calculate averages. avg_HS_bias = (double)sum_HS_bias/no_HS_lines; avg_MY_tput = (double)sum_MY_tput/no_MY_lines; avg_HS_tput = (double)sum_HS_tput/no_HS_lines; avg_MY_SNR = (double)sum_MY_SNR/no_MY_lines; avg_HS_SNR = (double)sum_HS_SNR/no_HS_lines; avg_MY_dfrq = (double)sum_MY_dfrq/no_MY_lines; avg_HS_dfrq = (double)sum_HS_dfrq/no_HS_lines; avg_MY_phs = (double)sum_MY_phs/no_MY_lines; avg_HS_phs = (double)sum_HS_phs/no_HS_lines; avg_MY_pwr = (double)sum_MY_pwr/no_MY_lines; avg_HS_pwr = (double)sum_HS_pwr/no_HS_lines; no_MY_num_ecc = no_MY_lines - no_MY_XX; //Calculate #s of numerical ecc percentages. no_HS_num_ecc = no_HS_lines - no_HS_XX; if((no_MY_num_ecc < 2) || (no_HS_num_ecc < 2)) {printf("Not enough data for statistics.\n"); return 0;} avg_MY_ecc = (double)sum_MY_ecc/no_MY_num_ecc; avg_HS_ecc = (double)sum_HS_ecc/no_HS_num_ecc; //Calculate standard deviations. sd_MY_bias = sqrt((double)sum_sq_MY_bias/(no_MY_lines - 1) - (double)no_MY_lines*avg_MY_bias*avg_MY_bias/(no_MY_lines - 1)); sd_HS_bias = sqrt((double)sum_sq_HS_bias/(no_HS_lines - 1) - (double)no_HS_lines*avg_HS_bias*avg_HS_bias/(no_HS_lines - 1)); sd_MY_tput = sqrt((double)sum_sq_MY_tput/(no_MY_lines - 1) - (double)no_MY_lines*avg_MY_tput*avg_MY_tput/(no_MY_lines - 1)); sd_HS_tput = sqrt((double)sum_sq_HS_tput/(no_HS_lines - 1) - (double)no_HS_lines*avg_HS_tput*avg_HS_tput/(no_HS_lines - 1)); sd_MY_SNR = sqrt((double)sum_sq_MY_SNR/(no_MY_lines - 1) - (double)no_MY_lines*avg_MY_SNR*avg_MY_SNR/(no_MY_lines - 1)); sd_HS_SNR = sqrt((double)sum_sq_HS_SNR/(no_HS_lines - 1) - (double)no_HS_lines*avg_HS_SNR*avg_HS_SNR/(no_HS_lines - 1)); sd_MY_dfrq = sqrt((double)sum_sq_MY_dfrq/(no_MY_lines - 1) - (double)no_MY_lines*avg_MY_dfrq*avg_MY_dfrq/(no_MY_lines - 1)); sd_HS_dfrq = sqrt((double)sum_sq_HS_dfrq/(no_HS_lines - 1) - (double)no_HS_lines*avg_HS_dfrq*avg_HS_dfrq/(no_HS_lines - 1)); sd_MY_phs = sqrt((double)sum_sq_MY_phs/(no_MY_lines - 1) - (double)no_MY_lines*avg_MY_phs*avg_MY_phs/(no_MY_lines - 1)); sd_HS_phs = sqrt((double)sum_sq_HS_phs/(no_HS_lines - 1) - (double)no_HS_lines*avg_HS_phs*avg_HS_phs/(no_HS_lines - 1)); sd_MY_ecc = sqrt((double)sum_sq_MY_ecc/(no_MY_num_ecc - 1) - (double)no_MY_num_ecc*avg_MY_ecc*avg_MY_ecc/(no_MY_num_ecc - 1)); sd_HS_ecc = sqrt((double)sum_sq_HS_ecc/(no_HS_num_ecc - 1) - (double)no_HS_num_ecc*avg_HS_ecc*avg_HS_ecc/(no_HS_num_ecc - 1)); sd_MY_pwr = sqrt((double)sum_sq_MY_pwr/(no_MY_lines - 1) - (double)no_MY_lines*avg_MY_pwr*avg_MY_pwr/(no_MY_lines - 1)); sd_HS_pwr = sqrt((double)sum_sq_HS_pwr/(no_HS_lines - 1) - (double)no_HS_lines*avg_HS_pwr*avg_HS_pwr/(no_HS_lines - 1)); printf("Number of reports on %s = %d\n",MY,no_MY_lines); //Print stats. printf("Number of reports on %s = %d\n\n",HS,no_HS_lines); printf("E(MY_bias) = %2.1f, sd(MY_bias) = %3.2f\n",avg_MY_bias,sd_MY_bias); printf("E(HS_bias) = %2.1f, sd(HS_bias) = %3.2f\n",avg_HS_bias,sd_HS_bias); printf("E(MY_tput) = %4.2f cps, sd(MY_tput) = %4.2f cps\n",avg_MY_tput,sd_MY_tput); printf("E(HS_tput) = %4.2f cps, sd(HS_tput) = %4.2f cps\n",avg_HS_tput,sd_HS_tput); printf("E(MY_SNR) = %4.2f dB, sd(MY_SNR) = %4.2f dB\n",avg_MY_SNR,sd_MY_SNR); printf("E(HS_SNR) = %4.2f dB, sd(HS_SNR) = %4.2f dB\n",avg_HS_SNR,sd_HS_SNR); printf("E(MY_dfrq) = %4.2f Hz, sd(MY_dfrq) = %4.2f Hz\n",avg_MY_dfrq,sd_MY_dfrq); printf("E(HS_dfrq) = %4.2f Hz, sd(HS_dfrq) = %4.2f Hz\n",avg_HS_dfrq,sd_HS_dfrq); printf("E(MY_phs) = %4.2f deg., sd(MY_phs) = %4.2f deg.\n",avg_MY_phs,sd_MY_phs); printf("E(HS_phs) = %4.2f deg., sd(HS_phs) = %4.2f deg.\n",avg_HS_phs,sd_HS_phs); printf("E(MY_ecc) = %4.2f%%, sd(MY_ecc) = %4.2f%%\n",avg_MY_ecc,sd_MY_ecc); printf("E(HS_ecc) = %4.2f%%, sd(HS_ecc) = %4.2f%%\n",avg_HS_ecc,sd_HS_ecc); printf("E(MY_pwr) = %4.2f%%, sd(MY_pwr) = %4.2f%%\n",avg_MY_pwr,sd_MY_pwr); printf("E(HS_pwr) = %4.2f%%, sd(HS_pwr) = %4.2f%%\n\n",avg_HS_pwr,sd_HS_pwr); printf("Error correcting capability exceeded %d time%c on MY link\n",no_MY_XX,(no_MY_XX==1)?' ':'s'); printf("Error correcting capability exceeded %d time%c on HS link\n\n",no_HS_XX,(no_HS_XX==1)?' ':'s'); printf("maximum \"MY\" throughput = %d cps, maximum \"MY\" SNR = %d dB\n",max_MY_tput,max_MY_SNR); printf("maximum \"HS\" throughput = %d cps, maximum \"HS\" SNR = %d dB\n\n",max_HS_tput,max_HS_SNR); printf("Modulation counts:\n"); printf("MY: %d BPSM %d QPSM %d 8PSM %d 8P2A %d 16P4A %d 16PSM\n",no_MY_BPSM,no_MY_QPSM,no_MY_8PSM,no_MY_8P2A,no_MY_16P4A,no_MY_16PSM); printf("HS: %d BPSM %d QPSM %d 8PSM %d 8P2A %d 16P4A %d 16PSM\n",no_HS_BPSM,no_HS_QPSM,no_HS_8PSM,no_HS_8P2A,no_HS_16P4A,no_HS_16PSM); fprintf(output1,"Number of reports on %s = %d\n",MY,no_MY_lines); fprintf(output1,"Number of reports on %s = %d\n\n",HS,no_HS_lines); fprintf(output1,"E(MY_bias) = %2.1f, sd(MY_bias) = %3.2f\n",avg_MY_bias,sd_MY_bias); fprintf(output1,"E(HS_bias) = %2.1f, sd(HS_bias) = %3.2f\n",avg_HS_bias,sd_HS_bias); fprintf(output1,"E(MY_tput) = %4.2f cps, sd(MY_tput) = %4.2f cps\n",avg_MY_tput,sd_MY_tput); fprintf(output1,"E(HS_tput) = %4.2f cps, sd(HS_tput) = %4.2f cps\n",avg_HS_tput,sd_HS_tput); fprintf(output1,"E(MY_SNR) = %4.2f dB, sd(MY_SNR) = %4.2f dB\n",avg_MY_SNR,sd_MY_SNR); fprintf(output1,"E(HS_SNR) = %4.2f dB, sd(HS_SNR) = %4.2f dB\n",avg_HS_SNR,sd_HS_SNR); fprintf(output1,"E(MY_dfrq) = %4.2f Hz, sd(MY_dfrq) = %4.2f Hz\n",avg_MY_dfrq,sd_MY_dfrq); fprintf(output1,"E(HS_dfrq) = %4.2f Hz, sd(HS_dfrq) = %4.2f Hz\n",avg_HS_dfrq,sd_HS_dfrq); fprintf(output1,"E(MY_phs) = %4.2f deg., sd(MY_phs) = %4.2f deg.\n",avg_MY_phs,sd_MY_phs); fprintf(output1,"E(HS_phs) = %4.2f deg., sd(HS_phs) = %4.2f deg.\n",avg_HS_phs,sd_HS_phs); fprintf(output1,"E(MY_ecc) = %4.2f%%, sd(MY_ecc) = %4.2f%%\n",avg_MY_ecc,sd_MY_ecc); fprintf(output1,"E(HS_ecc) = %4.2f%%, sd(HS_ecc) = %4.2f%%\n",avg_HS_ecc,sd_HS_ecc); fprintf(output1,"E(MY_pwr) = %4.2f%%, sd(MY_pwr) = %4.2f%%\n",avg_MY_pwr,sd_MY_pwr); fprintf(output1,"E(HS_pwr) = %4.2f%%, sd(HS_pwr) = %4.2f%%\n\n",avg_HS_pwr,sd_HS_pwr); fprintf(output1,"Error correcting capability exceeded %d time%c on MY link\n",no_MY_XX,(no_MY_XX==1)?' ':'s'); fprintf(output1,"Error correcting capability exceeded %d time%c on HS link\n\n",no_HS_XX,(no_HS_XX==1)?' ':'s'); fprintf(output1,"maximum \"MY\" throughput = %d cps, maximum \"MY\" SNR = %d dB\n",max_MY_tput,max_MY_SNR); fprintf(output1,"maximum \"HS\" throughput = %d cps, maximum \"HS\" SNR = %d dB\n\n",max_HS_tput,max_HS_SNR); fprintf(output1,"Modulation counts:\n"); fprintf(output1,"MY: %d BPSM %d QPSM %d 8PSM %d 8P2A %d 16P4A %d 16PSM\n",no_MY_BPSM,no_MY_QPSM,no_MY_8PSM,no_MY_8P2A,no_MY_16P4A,no_MY_16PSM); fprintf(output1,"HS: %d BPSM %d QPSM %d 8PSM %d 8P2A %d 16P4A %d 16PSM\n",no_HS_BPSM,no_HS_QPSM,no_HS_8PSM,no_HS_8P2A,no_HS_16P4A,no_HS_16PSM); //Write mysnrtpt setup file for gnuplot. fprintf(output4,"set pointsize 1\n"); fprintf(output4,"set rmargin 5\n"); fprintf(output4,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output4,"set xdata time\n"); fprintf(output4,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output4,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output4,"set title \"%s: E(tpt) = %4.2f cps E(SNR) = %4.2f dB\"\n",\ MY,avg_MY_tput,avg_MY_SNR); fprintf(output4,"plot '%s' using 1:4 title \"tpt\" with linespoints 1 1,\ '%s' using 1:5 title \"SNR\" with linespoints 3 3",fname2,fname2); //Write hssnrtpt setup file for gnuplot. fprintf(output5,"set pointsize 1\n"); fprintf(output5,"set rmargin 5\n"); fprintf(output5,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output5,"set xdata time\n"); fprintf(output5,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output5,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output5,"set title \"%s: E(tpt) = %4.2f cps E(SNR) = %4.2f dB\"\n",\ HS,avg_HS_tput,avg_HS_SNR); fprintf(output5,"plot '%s' using 1:4 title \"tpt\" with linespoints 1 1,\ '%s' using 1:5 title \"SNR\" with linespoints 3 3",fname3,fname3); //Write myphsecc setup file for gnuplot. fprintf(output6,"set pointsize 1\n"); fprintf(output6,"set rmargin 5\n"); fprintf(output6,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output6,"set xdata time\n"); fprintf(output6,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output6,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output6,"set title \"%s: E(phs) = %4.2f deg. E(ecc|<= 100) = %4.2f%%\"\n",\ MY,avg_MY_phs,avg_MY_ecc); fprintf(output6,"plot '%s' using 1:7 title \"phs\" with linespoints 1 1,\ '%s' using 1:8 title \"ecc\" with linespoints 3 3",fname2,fname2); //Write hsphsecc setup file for gnuplot. fprintf(output7,"set pointsize 1\n"); fprintf(output7,"set rmargin 5\n"); fprintf(output7,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output7,"set xdata time\n"); fprintf(output7,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output7,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output7,"set title \"%s: E(phs) = %4.2f deg. E(ecc|<= 100) = %4.2f%%\"\n",\ HS,avg_HS_phs,avg_HS_ecc); fprintf(output7,"plot '%s' using 1:7 title \"phs\" with linespoints 1 1,\ '%s' using 1:8 title \"ecc\" with linespoints 3 3",fname3,fname3); //Write mydfreq setup file for gnuplot. fprintf(output8,"set pointsize 1\n"); fprintf(output8,"set rmargin 5\n"); fprintf(output8,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output8,"set xdata time\n"); fprintf(output8,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output8,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output8,"set title \"%s: E(dfreq) = %4.2f deg.\"\n",MY,avg_MY_dfrq); fprintf(output8,"plot '%s' using 1:6 title \"dfreq\" with linespoints 1 1",fname2); //Write hsdfreq setup file for gnuplot. fprintf(output9,"set pointsize 1\n"); fprintf(output9,"set rmargin 5\n"); fprintf(output9,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output9,"set xdata time\n"); fprintf(output9,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output9,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output9,"set title \"%s: E(dfreq) = %4.2f deg.\"\n",HS,avg_HS_dfrq); fprintf(output9,"plot '%s' using 1:6 title \"dfreq\" with linespoints 1 1",fname3); //Write myhspwr setup file for gnuplot. fprintf(output10,"set pointsize 1\n"); fprintf(output10,"set rmargin 5\n"); fprintf(output10,"set timefmt \"%%H/%%M/%%S\"\n"); fprintf(output10,"set xdata time\n"); fprintf(output10,"set format x \"%%H/%%M/%%S\"\n"); fprintf(output10,"set xlabel \"[input: %s] time (HH/MM/SS) -->\"\n",fname0); fprintf(output10,"set title \"%s: E(pwr) = %4.2f%% %s: E(pwr) = %4.2f%%.\"\n",MY,avg_MY_pwr,HS,avg_HS_pwr); fprintf(output10,"plot '%s' using 1:9 title \"mypwr\" with linespoints 1 1,\ '%s' using 1:9 title \"hspwr\" with linespoints 3 3",fname2,fname3); return(0); } //End main(). /************************************************************************************* getline function: reads a line, returns its length. *************************************************************************************/ int getline(char *line, int max, FILE *iop) { if(fgets(line, max, iop) == NULL) return -1; else return strlen(line); } //End getline function. 
